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Abstract. We present the results of an extensive theo- 
retical investigation on the pulsation behavior of Bump 
Cepheids. We constructed several sequences of full ampli- 
tude, nonlinear, convective models by adopting a chemi- 
cal composition typical of Large Magellanic Cloud (LMC) 
Cepheids (Y=0.25, Z=0.008) and stellar masses ranging 
from M/Mq =6.55 to 7.45. We find that theoretical light 
and velocity curves reproduce the HP, and indeed close 
to the blue edge the bump is located along the descend- 
ing branch, toward longer periods it crosses at first the 
luminosity/ velocity maximum and then it appears along 
the rising branch. In particular, we find that the predicted 
period at the HP center is Php = 11.24 ± 0.46 d and that 
such a value is in very good agreement with the empiri- 
cal value estimated by adopting the Fourier parameters of 
LMC Cepheid light curves i.e. P H p = 11.2±0.8 d (Welch et 
al. 1997). Moreover, light and velocity amplitudes present 
a "double-peaked" distribution which is in good quali- 
tative agreement with observational evidence on Bump 
Cepheids. It turns out that both the skewness and the 
acuteness typically show a well-defined minimum at the 
HP center and the periods range from Php = 10.73 ±0.97 
d to Php = 11.29 ± 0.53 d which are in good agreement 
with empirical estimates. We also find that the models at 
the HP center are located within the resonance region but 
not on the 2:1 resonance line (P2/Po= 0.5), and indeed 
the P2/P0 ratios roughly range from 0.51 (cool models) to 
0.52 (hot models). 

Interestingly enough, the predicted Bump Cepheid 
masses, based on a Mass-Luminosity (ML) relation which 
neglects the convective core overshooting, are in good 
agreement with the empirical masses of Galactic Cepheids 
estimated by adopting the Baade-Wesselink method 
(Gieren 1989). As a matter of fact, the observed mass 
at the HP center -P 11.2 d- is 6.9 ± 0.9 M Q , while 
the predicted mass is 7.0 ± 0.45 M . Even by account- 
ing for the metallicity difference between Galactic and 
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LMC Cepheids, this result seems to settle down the long- 
standing problem of the Bump mass discrepancy. 

Finally, the dynamical behavior of a cool Bump 
Cepheid model provides a plain explanation of an ill- 
understood empirical evidence. In fact, it turns out that 
toward cooler effective temperatures the bump becomes 
the main maximum, while the true maximum is the bump 
which appears along the rising branch. This finding also 
supplies a plain explanation of the reason why the pul- 
sation amplitudes of Bump Cepheids present a "double- 
peaked" distribution. 

Key words: Stars: variables: Cepheids - hydrodynamics 
- Magellanic Clouds - stars: distances - stars: evolution - 
stars: oscillations 



1. INTRODUCTION 

More than seventy years ago Hertzsprung (1926) dis- 
covered that a subsamplc of Galactic classical Cepheids 
presents a relationship between the bump along the light 
curve and the pulsation period; the so-called "Hertzsprung 
Progression" was subsequently discovered among Mag- 
ellanic Clouds (MCs) and Andromeda Cepheids by 
Kukarkin & Parenago (1937), Shapley & Mckibben Nail 
(1940), and by Payne-Gaposchkin (1951,1954). The HP 
observational scenario was enriched by Joy (1937) and by 
Ledoux & Walraven (1958) who found a similar behavior 
in radial velocity curves. 

The empirical finger-print of the HP is the following: 
classical Cepheids in the period range 6 < P < 16 d show 
a bump along both the light and the velocity curves. This 
secondary feature appears on the descending branch of 
the light curve for Cepheids with periods up to 9 days, 
while it appears close to maximum light for 9 < P < 12 d 
and moves at earlier phases for longer periods. On the 
basis of this observational evidence this group of vari- 
ables was christened " Bump Cepheids" for avoiding to be 
mixed- up with "Beat Cepheids". In fact, the latter group 
refer to mixed-mode variables -i.e. objects in which two 
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or more modes are simultaneously excited- and therefore 
both the shape of the light curves and the pulsation ampli- 
tudes change from one cycle to the next, whereas Bump 
Ccpheids are single mode variables and their pulsation 
properties are characterized by a strong regularity over 
consecutive cycles. 

A more quantitative approach concerning Bump 
Ccpheids was originally suggested by Kukarkin & Pare- 
nago (1937), Payne-Gaposchkin (1947), and more recently 
by Simon & Lee (1981) who investigated the shape of the 
light curves by means of Fourier analysis. The last au- 
thors found that both the phase difference - fci- and the 
amplitude ratio -R21- show a sharp minimum close to the 
HP center. Following this approach several investigations 
have been already devoted to Fourier parameters of Galac- 
tic and Magellanic Cepheids. In particular, Moskalik et al. 
(1992, hereinafter MBM) suggested that the minimum in 
the Fourier parameters for Galactic Cepheids takes place 
at P H p = 10.0 ± 0.5 d, while Moskalik et al. (2000) by 
investigating a sample of more than 100 radial velocity 
curves found Php = 9.95 ± 0.05 d. At the same time, 
Welch et al. (1997) by investigating a large sample of 
Ccpheids in the LMC estimated that the minimum in the 
Fourier parameters is located at Php = 11.2 ±0.8 d. Thus 
supporting the shift of the HP center toward longer peri- 
ods originally suggested by Payne-Gaposchkin (1951) and 
strengthened by Andreasen & Petersen (1987) and by An- 
dreasen (1988). More recently Bcaulicu (1998) suggested 
that the HP center in LMC and in Small Magellanic Cloud 
(SMC) Cepheids is located at P H p = 10.5 ± 0.5 d and 
Php = 11.0±0.5d respectively. Since these three stellar 
systems are characterized by different mean metallicities, 
namely Z=0.02 (Galaxy), Z=0.008 (LMC), and Z=0.004 
(SMC), this empirical evidence seems to suggest that a de- 
crease in metallicity moves the HP center toward longer 
periods. 

Up to now, two distinct models have been proposed in 
the literature to explain the appearance of the HP among 
Bump Ccpheids; the echo model and the resonance model. 
The former was suggested by Whitney (1956) and dis- 
cussed by Christy (1968,1975) on the basis of Cepheid non- 
linear, radiative models. According to Christy, during each 
cycle close to the phases of minimum radius and before the 
phase of maximum expansion velocity a pressure excess is 
generated in the first He ionization region. This pressure 
excess causes a rapid expansion which in turn generates 
two pressure waves moving outward and inward. The lat- 
ter reaches the stellar core close to the phase of maximum 
radius, then reflects and reaches the surface one cycle later 
causing the appearance of the bump. The resonance model 
was suggested by Simon & Schmidt (1976, hereinafter SS) 
and is based on linear, adiabatic periods. In this theoreti- 
cal framework the bump would be caused by a resonance 
between the second overtone and the fundamental mode 
and it takes place when the period ratio between these 
two modes is close to 0.5. In particular, they suggested 



that the instability of the fundamental mode drives, due 
to a resonance, the second overtone instability. This ex- 
planation lies on the evidence that the nonlinear, radia- 
tive models constructed by Stobie (1969) show a bump 
along the radial velocity curves close to the resonance line 
P 2 /P = 0.5. 

Such an extensive observational and theoretical effort 
devoted to Bump Cepheids were not only aimed at un- 
derstanding the HP but also at providing independent es- 
timates of both the mass and the radius of these vari- 
ables. In fact, dating back to Christy (1968, 1975), Stobie 
(1969) and Fricke et al. (1972, hereinafter FSS) it was 
suggested that these two evolutionary parameters can be 
constrained on the basis of period and phase of the bump. 
A different method to estimate the mass, based on pe- 
riod ratios, was suggested by Petersen (1973). Note that 
mass determinations based on these two methods present 
a compelling feature: they are based on observables such 
as periods and phases of the bump which are not affected 
by systematic empirical uncertainties and therefore they 
are only limited by photometric accuracy. However, pul- 
sational masses based on these methods are, with few ex- 
ceptions (Carson & Stothers 1988, hereinafter CS), sys- 
tematically smaller than the evolutionary masses. This 
longstanding puzzle raised the so-called Bump mass dis- 
crepancy (see also Cox 1980) and at the same time sup- 
ported the use of a ML relation based on evolutionary 
models which include a mild or a strong convective core 
overshooting (Simon 1995; Wood 1998). Even though, the 
new radiative opacities settled down this vexata questio 
(MBM; Kanbur & Simon 1994; Simon & Kanbur 1994, 
hereinafter SK), recent linear (Buchler et al. 1996; Simon 
& Young 1997) and nonlinear (Wood et al. 1997) predic- 
tions for Cepheids in the MCs present a small discrepancy 
with the ML relations predicted by current evolutionary 
models. 

The main aim of this paper is to use up-to-date non- 
linear hydrodynamical models which include the coupling 
between pulsation and convection as well as canonical evo- 
lutionary masses, to account for the observed properties 
of Bump Ccpheids. In this investigation we focus our at- 
tention on the classical HP and refer to a forthcoming 
investigation (Bono et al. in preparation) a more detailed 
discussion on the physical mechanisms which trigger the 
appearance of the HP among Bump Cepheids and on 
other resonances recently proposed for both short (An- 
tonello & Poretti 1986; Buchler et al. 1996) and long- 
period Cepheids (Pel 1978; Antoncllo 1998). In order to 
supply a detailed theoretical scenario we constructed a 
fine grid of full amplitude, nonlinear, convective models by 
adopting stellar masses ranging from 6.4 to 7.6 M Q and 
a fixed chemical composition, namely Y=0.25, Z=0.008. 
We adopted this chemical composition, because accurate 
observational data on LMC Cepheids are currently avail- 
able in the literature. At the same time, this metal abun- 
dance can supply useful constraints on the intrinsic accu- 
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racy of our nonlinear models, and indeed recent nonlinear, 
Ccpheid models constructed by adopting similar treat- 
ments of the coupling between pulsation and convection 
show at intermediate metal contents very large pulsation 
dcstabilizations and it has been suggested that a power- 
ful dissipation mechanism was not properly included in 
current pulsation codes (Buchlcr 2000). 

In §2 we briefly recall the theoretical framework 
adopted for constructing Cepheid models and present non- 
linear observables predicted by these models as well as 
their light and velocity curves. The systematic behavior 
of both luminosity and velocity amplitudes inside the in- 
stability strip is investigated in §3. In this section we also 
present the skewness and the acuteness of light and veloc- 
ity curves and discuss the use of these parameters to mark 
the position of the HP center. New analytical relations 
which connect the period at the HP center to the stellar 
mass and to the effective temperature of Bump Cepheid 
models at minimum amplitude are also provided in this 
section. In §4 we discuss the location of these models in 
the HR diagram, and compare the HP center predicted 
by nonlinear models with linear and nonlinear resonance 
lines. The dynamical behavior of two Bump Cepheid mod- 
els located close to the HP center and to the red edge of 
the instability strip is investigated in detail in sect. 5. The 
main results of this investigation are summarized in §6 
together with a brief discussion on the observables which 
can supply tight constraints on the accuracy of theoretical 
models. 

2. Theoretical framework 

The theoretical framework adopted for constructing full 
amplitude, nonlinear, convective Cepheid models was al- 
ready described in Bono et al. (1998) and in paper I (Bono 
et al. 1999) and therefore it is not discussed further here. 
The reader interested in the physical assumptions adopted 
to account for the coupling between pulsation and convec- 
tion is referred to Buchlcr (2000). In paper III (Bono et al. 
in preparation) we found that the sequence of models con- 
structed by assuming M/M & =7.0 and Y=0.25, Z=0.008 
showed a bump along both light and velocity curves which 
appears at earlier pulsation phases when moving toward 
longer periods. Accordingly, the amplitudes of these mod- 
els showed a "double-peaked" distribution with the two 
maxima located close to the blue (hot) and to the red 
(cool) edge of the instability strip and a well-defined min- 
imum close to P « 10.7 d. On the basis of these features 
we concluded that this sequence of models undergoes the 
HP when moving from the blue to the red edge. A similar 
behavior in the distribution of radial velocity amplitudes 
was also found by CS in a nonlinear, radiative investiga- 
tion of Bump Cepheid models. However as noted by the 
same authors, the velocity amplitudes were systematically 
larger than the observed ones and the predicted periods 
at the HP center were 10% longer than expected. 



In order to supply a more detailed mapping of the HP 
inside the instability strip we implemented the sequence 
at M/Mq —7.0 with new series of models constructed by 
adopting a mass step of 0.15 Af©. Linear and nonlinear 
calculations were performed by adopting the same input 
physics (opacity, equation of state) and the same ML re- 
lation adopted in our previous investigations and based 
on evolutionary models which neglect the convective core 
overshooting (see paper I and III for more details). The 
new sequences were extended in mass until we found a 
well-defined minimum in the pulsation amplitudes inside 
the instability strip. For providing a robust relationship 
between the bump progression and the pulsation period 
we adopted a temperature step of 100 K throughout the 
instability strip and of only 50 K where the bump crosses 
the luminosity maximum and moves from the descending 
to the rising branch. We end up with six new series of 
models ranging in mass from 6.55 to 7.45 M© which show 
the typical behavior of the HP. After the initial pertur- 
bation of the static model (see paper I) the approach to 
the nonlinear limit cycle required a direct time integration 
ranging from 1500 to 5000 cycles. 

The input parameters of each sequence are listed in 
Table 1. This table gives the nonlinear pulsation prop- 
erties of the models which present a stable limit cycle. 
We typically adopted a temperature step of 100 K, there- 
fore the edges of the instability region can be estimated 
by increasing/decreasing the effective temperature of the 
hottest/coolest model by 50 K. The first three columns 
list for each model the stellar mass (solar units), the loga- 
rithmic luminosity (solar units), and the static effective 
temperature (K). Columns 4) and 5) give the nonlin- 
ear, fundamental period (d) and the logarithmic mean 
radius (solar units). The observables listed in the other 
columns are the following: 6) fractional radius oscillation, 
i.e. AR/R ph = (i? max - R min )/R ph where R ph is the pho- 
tospheric radius; 7) radial velocity amplitude (km s _1 ), 
i.e. Au = u max ~ u mm ; 8) bolometric amplitude (mag), 
i.e. AMboi = M^f - M£f; 9) logarithmic amplitude 
of static gravity, i.e. Alogg s = logg£ la * - logg™ m ; 10) 
logarithmic amplitude of effective gravity, i.e. Alogg c g — 
logg™£ x - logg™i n where g cS = GM/R 2 + du/dt; 11) tem- 
perature amplitude (K), i.e. AT = T max - T min where 
T is the temperature of the outer boundary; 12) effec- 
tive temperature amplitude (K), i.e. AT C = T c max - T c min 
where T e is derived from the surface luminosity, 13) total 
Kinetic energy (erg). The quantities listed in columns 4) 
to 11) refer, with the exception of the effective gravity, to 
the surface zone. 

Figs. 1-6, show the light and the velocity curves of the 
six new sequences for the labeled values of both mass and 
luminosities. Light and velocity curves of the sequence for 
M/M© =7.0 were already published in the ApJ on-line 
edition of paper III (Fig. llf). Note that light and veloc- 
ity curves plotted in these figures show neither spurious 
secondary features (bumps or dips) along the pulsation 
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cycle nor sudden jumps close to the phases of maximum 
compression. As already noted in paper I and III, this 
feature is a significative improvement in comparison with 
light and velocity curves predicted by nonlinear, radia- 
tive models (Christy 1975), and indeed theoretical curves 
were plotted without applying any running average or fil- 
tering process for smoothing the surface variations (Karp 
1975). This finding is even more relevant for investigating 
the HP. In fact, up to now nonlinear, radiative predic- 
tions were mainly based on radial velocity curves, since 
the light curves presented several spurious secondary fea- 
tures (CS; MBM). As a consequence, even though the bulk 
of observational data on the HP comes from light curves, 
theoretical insights were focused on the shape of radial 
velocity curves. 

Curves displayed in Figs. 1-6 show quite clearly that 
when moving from hotter to cooler effective tempera- 
tures the bump moves from the descending branch to the 
rising branch, i.e. an increase in the period moves the 
bump along the transition: descending branch — ► lumi- 
nosity maximum — > rising branch. Theoretical predictions 
in Figs. 1-6 disclose two key features: 

a) the bump along the light curves crosses the lumi- 
nosity maximum at shorter periods when compared with 
the velocity curves. This means that at fixed mass and for 
decreasing effective temperatures the Cepheid light curves 
start to show a flat-topped shape when the velocity curves 
still present the bump along the descending branch. This 
finding suggests that the lasting of the flat-topped phase in 
the luminosity curve is mainly governed by radius rather 
than by temperature variations. 

b) Along each sequence the model which attains the 
smallest velocity amplitude presents a flat-topped shape 
and the minimum in the luminosity amplitude coincides, 
within current temperature resolution, with the minimum 
in the velocity amplitude. On the basis of this evidence we 
will assume as period of the HP center -Pap- the period 
of the model which attains, along each sequence, the min- 
imum amplitude in both luminosity and velocity changes. 

3. HP Systematic behavior 

In order to provide a quantitative estimate of the change in 
the pulsation properties we take into account three differ- 
ent parameters, namely the amplitude, the acuteness, and 
the skewness of both light and velocity curves. Although 
the Fourier parameters supply more detailed information 
on the dynamical behavior along the pulsation cycle, we 
adopted these three parameters because we are mainly in- 
terested in the HP center and also because they can be 
safely estimated even if the light and the velocity curves 
are not perfectly sampled. 

According to Stellingwerf & Donohoe (1987) we define 
the acuteness as the ratio of the phase duration during 
which the magnitude is fainter than the median magni- 
tude -m me d = 0.5 * (m max + m m in)- to the phase du- 



ration during which it is brighter -</>t,- than m me( j, i.e. 
Ac = (l/0b) — 1. This quantity is a measure of the top- 
down asymmetry of the light curve and it decreases when 
the shape changes from sawtooth to flat-topped. At the 
same time, we define the skewness as the ratio of the phase 
duration of the descending branch to the phase duration 
of the rising branch -<p T -, i.e. Sk — (1/4> T ) — 1- This quan- 
tity is a measure of the asymmetry between rising and 
decreasing branch and it decreases when the slope of the 
rising branch becomes flatter. For symmetrical curves -not 
necessarily sinusoidal- both Sk and Ac attain values close 
to 1. 

Table 2 gives in the first two columns the mass and 
the effective temperature of each model, and in column 
3) the ratio of second overtone to fundamental period. 
Since the second overtones are pulsationally stable in this 
region of the instability strip, we adopted for this mode 
the linear, nonadiabatic period. Columns 4) to 6) list the 
amplitude, the acuteness and the skewness of visual light 
curves. Bolometric magnitudes were transformed into V 
magnitudes by adopting bolometric corrections by Castclli 
et al. (1997). Columns 7) and 8) give the acuteness and 
the skewness of radial velocity curves. 

Fig. 7 shows the luminosity amplitude in the V band 
-Ay- (top panel), the skewness (middle panel), and the 
acuteness (bottom panel) as a function of the logarithmic 
period. Sequences characterized by different stellar masses 
(see labeled values) are plotted using different symbols. 
Data plotted in Fig. 7 show that these three parame- 
ters typically present well-defined minima. We find that 
an increase in the stellar mass moves these minima to- 
ward longer periods and that the periods of these minima 
are located at Pap(Ay) = 11.24 ± 0.46 d, Pap{Ac) = 
11.17 ±0.48 d, and Pap(Sk) = 10.73 ± 0.97 d, where the 
errors give the standard deviations. The period of the HP 
center based on Sk minima is shorter than Php(^4v), and 
Pap(Ac) because the sequence for M/M Q =6.55 does not 
show a sharp minimum. 

These theoretical minima appear in very good agree- 
ment with the empirical determination based on the 
Fourier parameters <p2i and i? 2 i of a large sample of LMC 
Cepheids provided by Welch et al. (1997) i.e. Pap = 
11.2 ± 0.8 d and in reasonable agreement with the esti- 
mate, based on the same approach, provided by Beaulieu 
(1998) i.e. P H p = 10.5 ± 0.5 d. Here we note that the ob- 
served Php value was estimated as the period at which the 
Fourier parameters 4>2i and R21 of the light curves present 
a sudden jump, while the predicted Php values are the pe- 
riods at which Bump Cepheid models attain the minimum 
value in Ay, Sk, and Ac respectively. However, Stellingw- 
erf & Donohoe (1986) using adiabatic, one-zone pulsation 
models showed that the jump in the <f>2i parameter is cor- 
related with a local minimum in the Sk parameter. At 
the same time, Andreasen & Petersen (1987) in a detailed 
analysis of both amplitudes and Fourier parameters found 
that the jump in 4>2i and in R21 of LMC Bump Cepheids 
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Fig. 7. Luminosity amplitude in ffie^V band (top pariklj-? 
skewness (middle panel), and acuteness (bottom panel) as L 
a function of logarithmic period. Models characterized by 
different stellar masses (see labeled values) were plotted 
using different symbols. 



takes place at the same period, within observational un- 
certainties, at which the B photographic amplitudes, the 
skewness, and the acuteness attain their minimum value. 
As a consequence, we conclude that not only Pnp(Ay) 
but also Pnp(Sk) and Php(Ac) are robust indicators of 
the HP center. 

Data plotted in the top panel show that when moving 
from the blue to the red edge the luminosity amplitudes 
do not show a monotonic behavior with stellar mass. In 
fact, close to the blue edge an increase in the stellar mass 
causes a decrease in the pulsation amplitudes, whereas 
close to the red edge the amplitudes present an opposite 
behavior. Moreover, the luminosity amplitudes present a 



'objects 'ai'e ch'aiV'le'rized by 1 a'spreall rrrtrte 'puls'al'ion pe- 
riod and thereWe the masses and' the effective tempera- 
^S-CfPddrrved using the previous relations can be affected 
by larger uncertainties. In fact, by adopting the empirical 
value -P H p = 11.2 ±0.8 d- provided by Welch et al. (1997) 
we find M HP = 7.0 ± 1.0M Q , and T C HP = 5300 ± 600 K. 

Fig. 8 shows the same parameters of Fig. 7 but for the 
radial velocity curves. The behavior of AV, Sk, and Ac 
versus period are quite similar to the behaviors of the same 
parameters of the light curves. We find that the periods of 
the minima are now located at Pnp(AV) = 11.24 ± 0.46 
d, Pkp(Ac) = 11.29 ± 0.53 d, and P U p(Sk) = 11.27 ±0.49 
d respectively, and are once again in good agreement with 
empirical estimates. Although light and velocity curves 
show quite similar Sk values over the whole period range, 
the velocity curves present larger Sk values just before the 
HP center. This effect is due to the fact that the shape 
of the velocity curves becomes flat-topped at minimum 
amplitude whereas the light curves show the same shape 
at shorter periods (see §2). The comparison of predicted 
radial velocity amplitudes with observational data is dif- 
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ficult since spectroscopic measurement of LMC Cepheids 
are very scanty or do not cover the HP (Antonello 1998). 
However, radial velocity amplitudes for Galactic Cepheids 
(see Fig. 9 of paper III) collected by Cogan (1980) and by 
Bcrsier et al. (1994) show that the upper envelope presents 
a well-defined minimum across the HP center. The same 
outcome applies to the skewness, and indeed recent spec- 
troscopic measurements of Galactic Cepheids (Gorynya 
1998) show a sharp decrease close to logP « 1. This pro- 
vides a qualitative support to the predicted minima and 
to the "double-peaked" distributions shown in Fig. 8. 

Finally, we note that the P2/P0 ratios (see Table 2) 
of the models located at the HP center range from 0.517 
for the M/Mq =6.55 model to 0.512 for the M/M Q =7.45 
model. These period ratios were estimated by adopting 
the linear second overtone periods, therefore we can safely 



The SK^resWance lines were evaluated by adopting a 
^^wbonLrtaphson algorithm which for each given value of 
the effective temperature constructs a sequence of linear, 
nonadiabatic, radiative models by iterating on the values 
of mass and luminosity until both the resonance condition 
-P2/P0 =0.5- and the resonance center -P = 10.5, 11.0 d- 
are matched. A slightly different approach for estimating 
the resonance lines in a linear regime was suggested by 
Buchler et al. (1996) but their predictions were not in- 
cluded in Fig. 9, because the effective temperature values 
were not provided. 

Fig. 9 shows that the SK resonance lines are located 
at hotter effective temperatures when compared with the 
nonlinear HP centers and this discrepancy increases to- 
ward higher luminosities. On the other hand, our linear 
and nonlinear resonance lines are located at effective tem- 
peratures systematically cooler that the nonlinear HP cen- 
ters. The difference between our and SK resonance lines is 
not surprising, since our pulsation models are constructed 
by adopting different physical assumptions and chemical 
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Fig. 9. Distribution in the HR diagram of Bump Ccphcid 
models. Cepheid models with a bump along the descend- 
ing or the rising branch are plotted as triangles and 
squares. Cepheid models which attain along each sequence 
the minimum amplitude are plotted as filled circles. The 
dotted and the dashed lines refer to the linear P2/P0 =0.5 
resonance lines derived by SK assuming that the resonance 
center is at Po = 10.5 d and 11.0 d respectively. The long- 
dashed and the solid lines show our linear and nonlinear 
P2/P0 =0.50 resonance lines. 



composition. However, we note that for each given lumi- 
nosity the stellar masses adopted by SK are systematically 
smaller than our mass values and this difference typically 
increases when moving from lower to higher luminosities. 

The same difference appears in mass estimates, based 
on a similar approach, provided by Beaulieu & Sasselov 
(1997, and references therein), and indeed for LMC Bump 
Ccpheids at the resonance center P = 10.5 d they give a 



In this context it is noteworthy that stellar masses 
based on linear period ratios are, as already noted by 
MBM, SK, and by Buchler et al. (1996), systematically 
smaller than the observed ones. On the other hand, our 
nonlinear, convective models based on a canonical ML 
relation agree with empirical Galactic Ccphcid masses 
estimated by Gieren (1989) using the Baade-Wesselink 
method. In fact, by using the Gieren's Period-Mass re- 
lation (his relation 2) for Php = 11 -2 d we obtain an em- 
pirical Cepheid mass of M/Mq — 6.9 ± 0.9. This mass 
value is in satisfactory agreement, within the observa- 
tional uncertainties, with nonlinear predictions, and in- 
deed at the same period the theoretical masses range from 
6.55 to 7.45 Mq. We know that we are comparing theo- 
retical predictions for LMC Cepheids with observational 
data for Galactic Cepheids. However, current uncertain- 
ties on mass determinations are probably larger than the 
metallicity effect on Cepheid masses. In fact, preliminary 
theoretical results support the evidence that the Bump 
Cepheid masses at solar chemical composition are of the 
order of 6.5 ± 0.25 M (Ricci 1999). 
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Fig. 10. Light (top panel) and velocity (bottom panel) 
variations along two consecutive cycles. Solid and dashed 
lines refer to fundamental models constructed by adopting 
the same luminosity and effective temperature but differ- 
ent composition and mass value. 



As a result, predictions based on nonlinear, convec- 
tive models and on canonical evolutionary tracks settle 
down the long-standing conundrum of the Bump mass dis- 
crepancy. However a firm conclusion on this problem can 
be reached as soon as pulsation calculations constructed 
by adopting a noncanonical ML relation (i.e. based on 
evolutionary models which account for convective core 
overshooting) will become available. Note that accurate 
and homogeneous mass determinations based on the near- 
infrared surface brightness technique (Gieren et al. 1997) 
and on the infrared flux method (Fernley at al. 1989) can 
supply tight constraints on the Cepheid ML relation. Inde- 
pendent mass estimates can also be derived from orbits of 



tal model located close thShc HP center (P2/P0 =0.516, 
PftASfe^ 00 K ) of tnc sequence at M/M =7.0. Regions in 
which the radial velocity attains positive/negative values 
are plotted as solid/dotted lines. For avoiding mislead- 
ing interpretations of the amplitude variations the radial 
velocities were plotted, according to Aikawa & Whitney 
(1985), without applying any artificial shift or enhance- 
ment. 

Data plotted in Fig. 11 show that the outermost layers 
contract on a shorter time scale than the regions located 
just below them. In fact, the external layers reach the 
minimum contraction velocity at <j) = 0.23 and then at 
(j) = 0.40 -the phase of minimum radius- start to expand, 
whereas the underlying regions are still contracting. We 
measure phases with respect to the beginning of the sur- 
face contraction at maximum radius i.e. (j>(R = R max ) = 0. 
Soon after this phase the outermost layers undergo at first 
a rapid expansion reaching the maximum expansion veloc- 
ity at (j) = 0.54 and then a decrease in the radial velocity 
due to gravity. During this slow down phase -<f> w 0.65- 
the outermost layers interact with deeper layers (see the 
arrow) that are rapidly moving outward and reach their 
maximum expansion velocity at <f> « 0.74. The interaction 
between these two different dynamical behaviors causes 



Bono et al.: The Hertzsprung progression 



9 



an increase in the acceleration both toward the surface 
and toward the center. As a consequence in the outermost 
layers the radial velocity increases once again and then 
slows down showing a secondary maximum -the bump- at 
(j> = 0.85. 

The innermost layers start to expand around <f> — 0.45 
but their expansion is not in phase with the outermost lay- 
ers. In fact, the latter ones reach the maximum expansion 
velocity at phases during which the bulk of the envelope is 
still accelerating. This phase difference causes the forma- 
tion of a shock between the slowing down motion of outer- 
most layers and the outgoing expansion of deeper layers. 
As a consequence, the shocked region undergoes a strong 
compression which in turn generates pressure waves that 
propagate both toward the surface and toward the center 
of the stellar model. It is the outgoing pressure wave that 
causes, as already noted, the appearance of the bump, 
whereas the incoming pressure wave -the Christy wave- 
causes a short contraction phase in the innermost regions 
lasting from <fi = 0.7 to <j> = 0.85. Soon after the bounce 
at the stellar core at <fi = 0.95, the Christy wave at first 
moves rapidly outward and shortly delays the beginning of 
the contraction phase in the innermost regions and then 
limits the inward excursion of the layers crossed during 
its propagation out to the surface. This effect causes close 
to the phases of minimum radius a bounce between the 
overlying layers -that are rapidly contracting- and these 
layers that are contracting more slowly. The shock formed 
by this bounce generates two pressure waves: the outgoing 
wave delays the outward excursion of the region located 
just below the outermost layers, while the incoming wave 
triggers the Christy wave. As a consequence, it is the prop- 
agation out to surface of the Christy wave generated in 
the previous cycle that eventually causes the phase shift 
between the outermost layers and of the underlying re- 
gions, and in turn the appearance of the bump. 

A phase shift between the outermost layers and the 
bulk of the envelope was also noted by Karp (1975). How- 
ever, he found that the Christy wave approaches the sur- 
face close to velocity maximum and not close to the bump 
phase, and therefore the echo model could not completely 
explain the HP. Our results support the timing suggested 
by Karp but it seems that it is the Christy wave which 
causes the phase shift, and in turn the bump. However, be- 
fore a firm conclusion on the plausibility of this mechanism 
can be reached a detailed comparison between theoretical 
predictions and empirical data has to be provided. 

We now briefly discuss the dynamical behavior of a 
model located close to the red edge (P2/P0 =0.502, T c = 
5100 K) of the sequence at M/Mq =7.0, since it presents 
an interesting feature worth being investigated. Fig. 12 
shows the radial velocity as viewed from the stellar sur- 
face. The overall behavior is similar to the model located 
close to the HP center, but in this red model the formation 
and the incoming propagation of the Christy wave can 
be more easily identified. In fact, close to 4> ~ 0-5 an d 



log < R > /10 12 w 3.5 its appearance is marked by a nar- 
row region in which the velocity attains negative values 
and then moves inward reaching the base of the envelope 
at 4> ~ 0.8 (see the arrows). In this model the contraction 
of the overlying layers is more rapid than for the model 
close to the HP center, and indeed the Christy wave stops 
the expansion of the underlying layers. This fact and the 
evidence that the outermost layers start to expand at later 
phases causes the appearance along the surface velocity 
curve of a shoulder close to 4> = 0.5 -j- 0.55 (see the ar- 
row at the top of the envelope). Soon after these layers 
experience a rapid outgoing acceleration and this motion 
is almost in phase with the motion of the envelope and 
therefore the outermost layers do not interfere with the un- 
derlying regions since they are almost co-moving. It turns 
out that in this model the bump, i.e. what wc defined as 
the secondary maximum, is the main maximum while the 
true maximum is the bump which appears on the rising 
branch. 

This finding provides a straightforward explanation to 
what has been defined by FSS as a happy but also ill- 
understood circumstance. In fact, these authors in order to 
derive a linear relationship between period, phase of the 
bump and Bump Cepheid radii noted that from an obser- 
vational point of view it could not be firmly stated whether 
the primary or the secondary maximum was generated by 
the Christy wave. To overcome this thorny problem they 
suggested to measure the phase of the secondary maxi- 
mum if the bump is located along the descending branch 
and the phase of the primary maximum if the bump ap- 
pears along the rising branch. The dynamical behavior of 
the Bump Cepheid model located close to the red edge 
supplies a plain theoretical support to this far-sighted ob- 
servational choice. At the same time, current models also 
explain why light and velocity amplitudes show a "double- 
peaked" distribution when moving from the blue to the 
red edge of the instability strip. The maximum located 
close to the blue edge is typical of relatively short-period 
Cepheids (see paper III) , whereas the second one is caused 
by the HP. 

6. Summary and conclusions 

This paper presents the results of an extensive theoreti- 
cal investigation on Bump Cepheids. In order to provide a 
detailed analysis of the pulsation behavior of these ob- 
jects we computed several sequences of full amplitude, 
nonlinear, convective models at fixed chemical composi- 
tion (Y=0.25, Z=0.008). The models were constructed by 
adopting a mass step of O.15M and a temperature step 
of 100 K close to the instability edges and of 50 K close 
to the HP center. The main outcomes of this analysis are 
the following: 

1) theoretical light and velocity curves account for the 
HP, and indeed close to the blue edge the bump is located 
along the descending branch, at longer periods it crosses at 
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first the luminosity/velocity maximum and then appears 
along the rising branch. 

2) In a very narrow period range both light and veloc- 
ity curves show a flat-topped shape and their amplitudes 
attain a well-defined minimum. The predicted period of 
the minima -Php = 11-24 ± 0.46 d- is in very good agree- 
ment with the empirical value found by Welch et al. (1997) 
for LMC Cepheids i.e. P H p = 11.2 ±0.8 d. Moreover, light 
and velocity amplitudes present a " double-peaked" distri- 
bution which agrees with observational data. 

3) Both the skewness and the acuteness of light and ve- 
locity curves show a well-defined minimum at the HP cen- 
ter and their distributions are in good qualitative agree- 
ment with empirical estimates provided by Andreasen & 
Petersen (1987). The periods of the HP center are, within 
the uncertainties, in good agreement with observational 
values. 

4) The models at the HP center are located within 
the resonance region but not on the 2:1 resonance line 
(P2/P0 =0-5), and indeed the P2/P0 ratios range from 0.51 
to 0.52. 

5) Predicted Bump Cepheid masses, based on a ML 
relation which neglects the convective core overshooting, 
are in good agreement with the empirical masses of Galac- 
tic Cepheids estimated by adopting the Baade-Wesselink 
method (Gieren 1989). In fact, the observed mass range at 
the HP center -P = 11.2 d- is 6.9±0.9 M Q , while the theo- 
retical one is 7.0 ±0.45 M Q . Even if Galactic Ccpheids arc 
more metal-rich than LMC Cepheids, this result seems to 
settle down the long-standing problem of the Bump mass 
discrepancy. 

The results presented in this paper were mainly aimed 
at testing whether current hydrodynamical models which 
include the coupling between pulsation and convection to- 
gether with a ML relation based on canonical evolutionary 
models account for the HP in LMC Cepheids. It is clear 
that the theoretical scenario we developed seems to pro- 
vide a reliable description of several empirical facts. More- 
over, we confirm that Bump Cepheid models can supply 
fundamental constraints on stellar masses, and in turn on 
the ML relation of intermediate-mass stars. We also note 
that current Cepheid metal-intermediate, convective mod- 
els do show plausible pulsation properties. 

Two crucial topics need to be properly addressed be- 
fore firm conclusions on the HP can be reached. To vali- 
date present nonlinear Cepheid models a thorough theo- 
retical analysis of the Fourier parameters should be pro- 
vided together with a detailed comparison with observa- 
tional data available in the literature (Bono et al. in prepa- 
ration). At the same time, the dependence of the HP be- 
havior both on the chemical composition should be inves- 
tigated as well. New and accurate CCD data on light and 
velocity amplitudes as well as on the skewness and the 
acuteness of these curves can certainly improve the loca- 
tion of the HP center and feed the future investigations 
with robust observational constraints. 
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Fig. 1. Bolometric lightpewses (left panel) and radial velocity (right panel) curves alongiAsB consecutive cycles as a 
function of the pulsation phase. The curves plotted in this figure refer to a selection of models constructed at fixed 
mass and luminosity (see labeled values). The nonlinear periods (d), and the effective temperatures (K) are listed in 
left and right panels. Positive and negative values denote expansion and contraction phases, respectively. 
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Similarto Fig. ]f?ii«£Efor M/M ^=6.85 models. 
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Similar to Fig. ]F>hbu±Efor a selection of models at M/M Q =7.15 
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Similar to Fig. l=hba±Efor a selection of models at M/M Q =7.3. 
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Similar to Fig. ]F>nbu±Efor a selection of models at M/M Q =7.45 
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Fig. 11. Radial velocity variations as a function of mean radius and phase for the model located close to the HP center, 
surface at bottom. Solid and dotted line refer to positive and negative velocities. The arrow marks the interaction 
between inner and outer layers. 




Fig. 12. Radial velocity variations as a function of mean radius and phase for a model located close to the red edge, 
surface at bottom. The two arrows at the base of the envelope mark the approach and the ensuing reflection of the 
Christy wave at the stellar core, while the arrow at the surface marks the phases at which appears the shoulder. 
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Table 1: Nonlinear Observables 
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a Stellar mass (solar units). b Logarithmic luminosity (solar units). c Effective temperature (K). d Fundamental 
period (day). e Logarithmic mean radius (solar units), f Fractional radius variation. 9 Radial velocity amplitude 
(Kms -1 ). h Bolometric amplitude (mag). 1 Amplitude of logarithmic static gravity. J Amplitude of logarithmic 
effective gravity. k Surface temperature variation (K). 1 Effective temperature variation (K). m Total kinetic energy 
(erg). 



Tabic 2: HP Systematic Behavior 
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Table 2: Continued 
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a Stellar mass (solar units). 6 Effective temperature (K). c Second overtone 
to fundamental period ratio. d Visual amplitude (mag). e Acuteness 
of visual light curve. ' Skewness of visual light curve. 9 Acuteness of 
radial velocity curve. h Skewness of radial velocity curve. 1 The Ac and 
Sk values based on light curves were estimated by adopting the magnitude 
scale -the maximum is the minimum value-, while the Ac and Sk values 
based on velocity curves were estimated by adopting the theoretical notation 
i.e. the maximum value is attained during the outward excursion (positive 
values) . Both for light and velocity curves we assumed that the phase is zero 
at maximum. 



Tabic 3: Input parameters and nonlinear results. 
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0.25 
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0.008 
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Au 


37.8 


38.8 
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0.71 


0.65 



a Y and Z are the helium and the metal content 
respectively. The other symbols and their units 
arc the same as in Tables land 2. b SK model. 
c Paper III model. 



